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ABSTRACT 



We present 12 CO (1-0) observations in the central 4.5 kpc (1') of the Ha/Radio 
lobe galaxy NGC 3079 with the Nobeyama Millimeter Array. The molecular gas 
shows four components: a main disk, spiral arms, a nuclear disk, and a nuclear 
core. The main disk extends along the galaxy major axis. We detected its central 
2 kpc radius, while its full extent is beyond our spatial coverage. Molecular gas 
is smoothly distributed in the main disk, having a gas mass of 5 x 1O 9 M within 
the central ~ 2 kpc radius. The spiral arms are superimposed on the main disk. 
Abrupt velocity changes of up to ~ 200 km s -1 are observed along the spiral arms 
in S-shaped twists of isovelocity contours and double velocity-peaked features on 
the spectra. The nuclear disk with ~ 600 pc radius appears in position- velocity 
(PV) diagrams, having an intense concentration of molecular gas. Its appearance 
on PV diagrams is indicative of oval motions of the gas, rather than circular. 
The nuclear disk and spiral arms form the so-called "figure-of-eight" pattern on 
a PV diagram. The nuclear core is more compact than our current resolution 
(2" = 150 pc), and has a gas mass of 3 x 1O 8 M within the central 150 pc. Though 
it is unresolved, the nuclear core shows a very high velocity ~ 200 km s -1 even 
at the radius of ~ 100 pc on the PV diagram. 

We propose a model that NGC 3079 contains a weak bar. The weak bar 
model explains the observed features of the main disk, spiral arms, and nuclear 
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disk. The main disk and spiral arms result from gaseous Xi-orbits and associated 
crowding respectively. The nuclear disk arises from gaseous x 2 -orbits. The gas 
concentration in the nuclear disk could be explained by the expected gas-fueling 
mechanism: the gas on rri-orbits flows along spiral arms (or offset shocks), collid- 
ing with the gas on X2-orbits, and accumulating onto the nuclear disk. Assuming 
that the gas moves nearly along the spiral arms which run perpendicular to the 
line-of-sight, the pattern speed of the bar is estimated to be 55 ± 10 km s _1 kpc -1 . 
The high velocity of the nuclear core cannot be explained by our model for a bar. 
Thus we attribute it to a central massive core with a dynamical mass of 1O 9 M 
within the central 100 pc. This mass is three orders of magnitude more massive 
than that of a central black hole in this galaxy. 

Subject headings: galaxies: active — galaxies: ISM — galaxies: kinematics and 
dynamics — galaxies: spiral — galaxies: structure 

1. INTRODUCTION 

NGC 3079 is a nearly edge-on (i = 77°) SBc galaxy at a distance of 15.6 Mpc (1 arcsec 
corresponds to 76 pc; Sofue et al. 1999). 

Its nucleus is classified as LINER (Heckman 1980) or Seyfert 2 (Ford et al. 1986; Ho, 
Fillipenko & Sargent 1997; Sosa-Brito, Tacconi-Garman & Lehnert 2001) from the optical 
emission spectrum, and shows strong H2O maser emission (Henkel et al. 1984; Nakai et al. 
1995; Trotter et al. 1998). The overall velocity distribution of H 2 masers suggests the 
presence of a binding mass of ~ 1O 6 M at the center (Trotter et al. 1998), possibly a central 
supermassive black hole, which has been found in many AGNs (Miyoshi et al. 1995; Wandel, 
Peterson & Malkan 1999; Ishihara et al. 2001). Parsec-scale nuclear jets observed in the 
radio continuum may be outflow from the central compact object (Irwin & Seaquist 1988; 
Sawada-Satoh et al. 2000). 

Kiloparsec-scale significant outflow along the galaxy minor axis is observed as lobes in 
the radio continuum (Duric & Seaquist 1988), Ha emission (Ford et al. 1986; Veilleux et 
al. 1994; Cecil et al. 2001, see Figure 1), and X-ray emission (Fabbiano, Kim & Trinchieri 
1992; Pietcsh, Trinchieri & Vogler 1998). Optical spectroscopy shows gas motions with a 
velocity range of ~ 2000kms _1 across the lobes, and unusually high [NII]/Hct line ratios 
which indicate the presence of shocks (Veilleux et al. 1994). This type of large-scale outflow 
is often attributed to starburst activity (Heckman, Armus & Miley 1990). However there are 
some arguments against this starburst model for NGC 3079 (Hawarden et al. 1995). Duric 
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& Seaquist (1988) present an alternative model that a wind flow from the nucleus could be 
directed toward the galaxy minor axis by interaction with dense gas surrounding the nucleus. 
This interaction produces shocks which explain the observed strength of H 2 emission from 
lower vibration transitions by collisional excitation (Hawarden et al. 1995). 

At the root of the kiloparsec-scale outflow, there is a dense molecular disk with a radius 
of a few kiloparsec (Young, Claussen & Scoville 1988; Sofue & Irwin 1992; Irwin & Sofue 
1992; Sofue et al. 2001). Based on their interferometry observations of CO (J = 1—0) 
emission, Sofue & Irwin (1992) found an intense concentration of gas in the molecular disk, 
spiral arm features visible on a position-velocity diagram, and a central compact core which 
was unresolved at their resolution (4"). Recently Sofue et al. (2001) have confirmed these 
features, and further resolved and classified the innermost region into the nuclear molecular 
disk and ultra-high-density core, based on their appearance on intensity maps and position- 
velocity diagrams. [In this paper, we will re-classify the features based on our new data 
(§3) and on theoretical considerations (§4).] The molecular disk was also observed in HCN 
and HCO + (Kohno et al. 2001), which measure the concentration of high density molecular 
clouds toward the nucleus. 

The disk of NGC 3079 is rotating with the north-side approaching and south-side reced- 
ing. A i^'-band image shows spiral arms in the disk, forming an inverted S'-shaped pattern 
on the sky (Veilleux, Bland-Hawthorn & Cecil 1999). If trailing spiral arms are assumed, 
the west of NGC 3079 is near-side, which is consistent with the entire dust lane morphol- 
ogy. NGC 3079 is a far-infrared (FIR) luminous galaxy listed in the IRAS bright galaxy 
catalog (Soifer et al. 1989). An ISO 90/im map shows that most of the FIR emission is 
produced by dust heated by stars in the entire galaxy disk, rather than the nucleus (Perez 
Garcia, Rodriguez Espinosa & Fuensalida 2000). The HI disk extends more widely than the 
optical disk and shows warps at the outskirts (Irwin & Seaquist 1991). These warps may 
originate from an interaction with a nearby companion NGC 3073. The HI velocity field, 
however, shows a regular rotating-disk pattern. The radial emission profile is well fitted by 
an exponential function, while it shows a sharp drop in the central region (< 50"; Irwin & 
Seaquist 1991) which is produced by absorption against strong continuum emission (Duric 
& Seaquist 1988). The CO disk is embedded in this region of the absorbed HI (Young, 
Claussen & Scoville 1988; Sofue & Irwin 1992; Sofue et al. 1999), and coincident with a void 
of Ha line emission (Cecil et al. 2001). 

Though NGC 3079 is classified as a barred galaxy (SBc), the stellar bar is hardly 
confirmed in optical/infrared photographs in this nearly edge-on galaxy. There is however 
some evidence for the presence of a bar. NGC 3079 has a stellar bulge which takes a so- 
called peanut-shape (Shaw, Wilkinson & Carter 1993), whose possible origin is a vertical 
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instability of rotating disk-stars in a bar potential (Combes et al. 1990). Veilleux, Bland- 
Hawthorn & Cecil (1999) fitted an oval orbit model to their Ha velocity field, and concluded 
that bar streaming motions with moderately eccentric orbits (e = b/a ~ 0.7) aligned along 
P. A. = 130° intrinsic to the disk (P. A. = 163° on the sky 5 ) are satisfactory to match the 
observations. 

In this paper, we present our recent 12 CO (J = 1 — 0) observations of the central 1' 
of NGC 3079 using the Nobeyama Millimeter Array (NMA). We describe the observations 
and data analysis in §2. A part of the data has been published in Sofue et al. (2001), but 
this paper presents our new analyses including new observations. The main features in the 
molecular disk of NGC 3079 are presented in §3. Gas dynamics in the molecular disk are 
discussed in §4. A weak bar and central massive core explain all the observed features. The 
central rotation curve and mass are derived in §5. We summarize our conclusions in §6. 



2. OBSERVATIONAL DATA 

2.1. NMA 12 CO (J = 1 - 0) Observations 

Our aperture synthesis observations of the 12 CO ( J = 1 — 0) emission from NGC 3079 
were obtained with the Nobeyama Millimeter Array (NMA) at the Nobeyama Radio Obser- 
vatory (NRO) 6 , between 2000 January and 2001 April for a single pointing center at (0:1950, 
^1950) — (9 h 58 m 35 s .02, +55°55'15".40). We made the observations with three available con- 
figurations (AB, C, and D); when combined, the visibility data cover projected baselines 
from 10 to 351 m. The NMA consists of six 10 m antennas, providing a FWHP of about 65" 
at 115 GHz. The antenna size limits the minimum projected baseline, restricting the largest 
detectable size of objects to about 54". Tunerless SIS receivers have receiver noise temper- 
atures of about 30 K in double sideband, and typical system noise temperatures of about 
400 K in single sideband. Digital spectro-correlators (Okumura et al. 2000) have two spec- 
troscopic modes; we used a mode covering 512 MHz (1331kms _1 ) with 2 MHz (5.2kms _1 ) 
resolution. We observed the quasar 0954+556 every 20 minutes for gain calibration, and 
3C279 (or 3C273) for bandpass calibration. Absolute flux scales (0.68 Jy at 115 GHz for 
0954+556, uncertain to ~ 20%) were measured three times in 2000 and once in 2001; no 



5 The bar position angle on the sky (P.A. s k y ) is calculated from that intrinsic to the disk (P.A.; n t), by 
"tan(P.A. ga i — P.A. s ky) = tan(P.A. ga i — P.A.i nt ) cos i" using their adopted position angle and inclination of 
the disk. 

6 Nobeyama Radio Observatory is a branch of the National Astronomical Observatory, operated by the 
Ministry of Education, Culture, Sports, Science and Technology, Japan. 
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significant flux variation was observed. 

The raw visibility data were calibrated for complex gain and passband with the UVPROC- 
II package developed at NRO, and mapped with the NRAO/AIPS package. We applied the 
CLEAN procedure with natural weighting for each velocity channel, and obtained a three- 
dimensional (RA, DEC, Vlsr) data cube. Parameters of the cube are listed in Table 2. 

Figure 2 displays zeroth- and first-moment maps in the central 20" x48" region (1.5 kpc x 
3.7 kpc) of NGC 3079, while Figure 3 shows channel maps with an interval of 10.4 km s -1 
in the same region. The significant emission (> 3a; la = 12mJybeam _1 ) is detected 
in 63 adjacent channels within the velocity range of Vlsr = 821 — 1467 km s -1 (AV = 
646 km s" 1 ). No primary beam correction has been applied in these maps. Figure 4 displays a 
position- velocity diagram (PV diagram) along the optically defined major axis (P.A.= 165°), 
integrated along the minor axis. The almost entire CO emission in Figure 2 (left) falls in 
the slit-width (12"). The axes are labeled relative to the dynamical center, and systemic 
recession velocity of the galaxy (derived in § 3.1.2; see Table 3). 

Compared with single-dish observations from IRAM 30 m and FCRAO 14 m telescopes 
(Braine et al. 1993; Young et al. 1995), our NMA cube recovers about 87% of the total CO 
line flux in the central 23" and 67% in 45", which means that the CO emission is concentrated 
in the central region. We made no correction for the missing flux in the following discussion. 
Figure 5 compares a spectrum from our NMA cube with that from the IRAM observations 
(Braine et al. 1993). In order to obtain the fluxes and spectrum comparable to the single-dish 
ones, the NMA cube was corrected for the primary beam response of antennas, convolved 
with the single-dish beam (23" for IRAM and 45" for FCRAO), and sampled at the pointing 
center of the single-dish observations. 

2.2. Supplied HST Data 

We obtained images from the Hubble Space Telescope (HST) archive (P.I. G. Cecil; see 
Cecil et al. 2001). The images were taken with WFPC2 and two filters, F814W and F658N, 
which correspond to the /-band and Ha + [N II] filters in ground telescopes. The central part 
of NGC 3079 lies on a WFC chip rather than the PC chip. Cosmic-ray hits were removed and 
images are combined with the IRAF/STSDAS package. The absolute position of the HST 
images is calibrated using the USNO-A2.0 catalog (Zacharias et al. 2000), and is accurate 
to about 0".5. The derived images are presented in the right-hand panel of Figure 1. 
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3. RESULTS 

Our maps and PV diagram show four distinct components: a main disk, spiral arms (or 
offset ridges), a nuclear disk, and a nuclear core. Figure 6 shows a schematic illustration of 
these four components on a PV diagram. We describe each of the four in order of decreasing 
radius in §3.1-3.4, and compare them with previous results in §3.5. 

3.1. The Main Disk 

3.1.1. Smooth Gas Distribution on the Main Disk 

The zeroth-moment map of Figure 2 displays a disk with an extent of 45" x 15" (3.4 kpc x 
1.1 kpc), elongated along the optically-defined major axis of the galaxy (P.A. = 165°). Since 
our synthesis observations have an intrinsic maximum limit of detectable scale (54"), larger 
components detected in single-dish observations 7 are not covered in this map, and result 
in the missing flux (our total flux recovery is about 67% in 45"). The inclination of the 
galaxy (77°) is large, but enough to show the emission distribution on the disk, because 
the vertical sizes of molecular disks are usually thin (150 pc for the Galaxy from Clemens, 
Sanders & Scoville 1988; 230 pc for NGC 891 from Handa et al. 1992) in comparison with the 
full extents of molecular disks (3.4 kpc or 45" for NGC 3079 from the zeroth-moment map). 
The CO emission is smoothly distributed on the full disk, showing no emission deficient 
regions often found in molecular disks (see Sakamoto et al. 1999a; Regan et al. 2001). The 
emission distribution is nearly axisymmetric, except for slight enhancements owing to spiral 
arms (see § 3.2). The main disk coincides with the void of HII regions (Cecil et al. 2001, see 
also Figure 7 left), and with the region where the HI gas is observed in absorption (Irwin 
& Seaquist 1991) against strong radio continuum emission (Duric & Seaquist 1988). The 
center of the main disk coincides with the root of the Ha-lobe. 



3.1.2. Kinematics of The Main Disk 

The first-moment map (Figure 2 right) shows a rotating disk with the northside ap- 
proaching and southside receding. The S'-shaped twists of isovelocity contours should be 
attributed to spiral arms (§ 3.2). The entire velocity field is almost perfectly bisymmetric, 



7 Recently, NRO 45 m telescope observations found that the full extent of CO emission in NGC 3079 is 
about 80" in radius (Yamauchi, private communication). 
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indicating regular motions in the molecular disk. 

We obtained kinematical parameters from the first-moment map using the AIPS / GAL 
package. The dynamical center, position angle and inclination are first determined using 
the Brandt rotation curve model (Brandt 1965). Then we fixed the above parameters, fit- 
ted tilted ring models with constant velocities to the first-moment map, and obtained the 
systemic recession velocity (LSR velocity). The results are listed in Table 3. The dynam- 
ical center coincides spatially with the emission centroid: (0:1950, ^1950) — (09 h 58 m 35 s .00, 
+55 d 55 m 15 s .80). The position angle and inclination are in good agreement with the results 
from optical isophotes (Sandage & Tamman 1981; Veilleux, Bland-Hawthorn & Cecil 1999) 
and from kinematics of the HI gas (Irwin & Seaquist 1991). The systemic recession velocity 
is consistent with estimates from optical spectroscopic data (e.g. 1150 km s' 1 from Veilleux, 
Bland-Hawthorn & Cecil 1999), but slightly exceeds the value from HI data (1124 km s" 1 
from Irwin & Seaquist 1991). This difference might come from a warp of the HI disk (Irwin 
& Seaquist 1991). We adopt the values from large-scale isophotes for position angle and 
inclination, and the derived values for galaxy center and systemic recession velocity in the 
rest of our analysis. These choices do not affect the following discussion. 

The PV diagram (Figure 4) shows the main disk as two ridges, symmetrically extending 
southeast and northwest from the velocity peaks around the center (see also Figure 6). The 
near-perfect symmetry in the diagram again indicates regular motions of the molecular gas 
in the disk. The velocity of 220 km s -1 at the radius of 26" (2.0 kpc) indicates the dynamical 
mass of Md yn = 2.2 x 1O 1O M . The disk is also confirmed in the channel maps (Figure 3); 
the approaching (NW) side appears in the channels of Vlsr = 821 — 1029 km s -1 , while the 
receding (SE) side is found in Vlsr = 1258 — 1467 km s -1 . 



3.1.3. Mass and Surface Densities of the Molecular Gas 

Masses of molecular gas Mh 2 are estimated from total CO-line flux S'co; CO-to-H 2 
conversion factor X co , and galaxy distance D by 

~ 72xl0 3 f — X( Sc ° ] 

x ( ^£2 

V 1.8 x 10 20 cm- 2 [Kkms- 1 ]- 1 

Assuming the hydrogen mass fraction (0.707 from Cox et al. 2000), the total gas mass, 
including He and other elements, becomes M gas = 1.41M H2 . Face-on surface densities of 
molecular gas are calculated from the integrated CO line intensity Iqq&V ', galaxy inclination 
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Jy km s -1 arcsec -2 

x ( ^£2 

Vl.8 x 10 20 cm- 2 [Kkms- 1 ]- 1 

The surface density of total gas also becomes E gas = 1.41E H2 . The surface density es- 
timation is independent of galaxy distance. We adopt the conversion factor of Xco = 
1.8 x 10 20 cm~ 2 [Kkms -1 ]" 1 from observations in the Galaxy (Dame, Hartmann and Thad- 
deus 2001), while X co could be smaller in galactic center regions (Arimoto, Sofue and 
Tsujimoto 1996). 

The total flux within the radius of 26" is S co = 1-9 x 10 3 Jykms" 1 , which corresponds 
to M gas = 4.6 x 1O 9 M for the galaxy distance D = 15.6 Mpc. The ratio of dynamical mass 
(2.2 x 1O 1O M from §3.1.2) to the total molecular gas mass M gas /M dyn is about 21% within 
the radius of 2kpc. The peak integrated intensity is Iqo^V = 1.3 x 10 2 Jybeam -1 kms -1 
(or 3.8 x lO^kms -1 ) at (a 1950 , 5 1950 ) = (09 h 58 m 35 s .00, +55°55'15 // .80); the corresponding 
face-on surface density becomes E gas = 2.5 x 1O 3 M pc -2 . No correction for missing flux is 
applied. 

Figure 8 shows the radial profile of the CO emission. Integrated intensities are az- 
imuthally averaged in each annulus (Ar = 0.5") with corrections for inclination (77°) and 
primary beam attenuation. The outskirts of the profile (r > 5") is well fitted by an exponen- 
tial (dotted line). At the innermost region, the CO emission exceeds this exponential curve, 
indicating a high CO concentration at the central region. Surface densities of molecular gas 
are labeled at the right-hand axis. The observed value of E H2 = 1000 M pc~ 2 at r ~ 5" 
(380 pc) is an order of magnitude higher than the typical value for the central regions of 
nearby galaxies (Sakamoto et al. 1999b). 

The HI gas surface density also follows an exponential profile E^ e~ r ^ hm where hm = 
118" (9.0 kpc) and = 8.8 x 10M o p C - 2 (Irwin & Seaquist 1991, from their Table 4). In 
comparison with this HI profile, our H 2 profile (dotted line) takes the same surface density 
at nearly the edge of our detected molecular disk (r ~ 23" or 1.7 kpc), rising more steeply 
toward the galaxy center (hco ~ 7.2" or 550 pc). This indicates that the gas in NGC 3079 
is highly concentrated in the central molecular disk. 



(— ) = 3.0xl0 2 cos*( 
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3.2. The Spiral Arms/Offset Ridges 

The zeroth-moment map (Figure 2 left) shows arm-like enhancements of emission, su- 
perimposed on the main disk; two ridges bisymmetrically run northwest and southeast from 
the galaxy center (most obvious within the central ~ 10" radius), forming an inverted S- 
shaped pattern. These features resemble (trailing) spiral arms, or offset ridges often found 
at the leading edge of stellar bars. [Note though these enhancements do arise from the 
spiral arms, they are slightly offset from the actual spiral arms (see Figure 7), because this 
projected map is contaminated by emission from other components.] The S-shaped twists 
of isovelocity contours on the first-moment map (Figure 2 left) coincide with the two arms, 
indicating non-circular motions around the spiral arms. Similar inverted S'-shaped spiral 
arms have been found in a i^'-band image (Veilleux, Bland-Hawthorn & Cecil 1999). 

Figure 9 shows line-profiles of CO emission at 5 x 7 points on a 3"-spacing grid (crosses 
on the left-hand panel), centered at the dynamical center (Table 3). Remarkable double- 
peaked features appear along the spiral arms: arrows point to examples of the double-peaks. 
The northwest arm corresponds to peaks with higher velocities, while the southeast arm cor- 
responds to peaks with lower velocities. The other peaks come from the main disk. These 
distinct velocities between the main disk and spiral arms (typically AV^ ~ 200 km s" 1 ) pro- 
duce the S'-shaped twists in the first-moment map. Similar double-peaked features have been 
noticed from Ha-spectroscopic data (Filippenko & Sargent 1992; Veilleux, Bland-Hawthorn 
& Cecil 1999), though their low velocity resolution did not allow them to distinguish between 
the main disk and spiral arms. 

The spiral arms are also confirmed on the PV diagram (Figure 4; see also Figure 6). Two 
ridges run northwest and southeast with lower velocities than those of the main disk. These 
two ridges and a central ridge of the nuclear disk (§3.3) show a so-called "figure-of-eight" 
pattern (or "tilted X" pattern). Similar patterns are often found in edge-on barred galaxies 
(Handa et al. 1990; Laine et al. 1990; Garcia-Burillo & Guelin 1995; Merrifield & Kuijken 
1999). In the velocity range of AV^ = ±150 kms -1 (V sys = 1147kms _1 ), the spiral arms are 
more dominant in emission than the main disk. 

This is also evident in the channel maps (Figure 3). The spiral arms are confirmed at 
V = 1010 — 1280 km s -1 ; emission from the northwest and southeast arms appear at the 
northwest and southeast edges of maps at 1010 and 1280 km s^ 1 , approaching the innermost 
core with increasing and decreasing velocities, respectively. Figure 7 (right-hand panel) 
shows contours of a zeroth-moment map made in the velocity range of AV^ = ±100 km s" 1 
centered at the systemic recession velocity (superimposed on the Ha + [Nil] map). The 
spiral arms are rooted to the nucleus, forming an inverted S-shaped pattern on the sky. 



-10- 



3.3. The Nuclear Disk 

The nuclear disk shows a sign of an existence on the PV diagram (Figure 4): remarkable 
velocity-peaks appear symmetrically at (—7", +250 km s -1 ) and (+7", —250 km s -1 ). Figure 
10 shows PV diagrams with 1 "-width and offsets of +2", 0" and —2" along the galaxy minor 
axis. The cuts are along the galaxy major axis. These diagrams reveal that the two peaks are 
on a single kinematical sequence: ridges are running and connecting the two peaks through 
three panels (see also Figure 6). This sequence seems to be rotating in this innermost region 
and forming a nuclear disk. The nuclear disk is obviously distinct from other velocity-peaks 
in the innermost region (at ±1" in Figure 4), because their corresponding peaks appear only 
in the diagram with 0"-offset; we thus think that they are consequences of the nuclear core 
(see § 3.4). The very high CO concentration within the nuclear disk (r < 7"; Figure 4 and 
8) indicates its distinct nature from the main disk. We conclude that the component which 
appears as the peaks and ridges in the PV diagrams is the nuclear disk. This component is 
rotating around the galaxy center with the northside approaching and southside receding. 

In the three panels of Figure 10, the ridges corresponding to the nuclear disk are not 
on a straight line connecting their two velocity peaks: the ridges are bent downward and 
upward in the +2" and —2" offset diagrams respectively, while the ridge twists upward (at 
the left-hand side) and downward (right-hand side) in the 0"-offset diagram. The deviations 
from the straight line amount up to ~ 50 km s _1 . These symmetrical bends and twist indicate 
that the gas in the nuclear disk takes oval orbits rather than circular (see § 4). 

The nuclear disk is also seen in the channel maps (Figure 3). The two velocity-peaks 
on the PV diagrams appears in the channels of V L sr = 822 — 926 km s^ 1 as an emission peak 
at around 7"-north from the center, and in Vlsr = 1384 — 1467 km s -1 at around 7"-south. 
In other channels, strong emission from the main disk and nuclear core contaminates the 
nuclear disk emission. 



3.4. The Nuclear Core 

The nuclear core is most clearly seen in the PV diagram (Figure 4), which appears as 
two velocity-peaks at (—1", +240 km s" 1 ) and (+1", —240 km s" 1 ). The distinct kinematics 
of this core from the other components is evident in a series of PV diagrams which have 
positional offsets along the galaxy minor axis (Figure 10). Among the three diagrams with 
+2", 0" and — 2"-offsets, the nuclear core appears only in the central one. The high velocities 
indicate a large dynamical mass in this nuclear region if the gas is pure-circularly rotating; 
we will discuss this in § 4.4. 
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The zeroth-moment map (Figure 2) shows the nuclear core as an intense CO concen- 
tration at the very center of the galaxy. It is also confirmed in channel maps (Figure 3) in 
almost all the velocities (Vlsr = 863 — 1436 km s -1 ). The central gas density of the nuclear 
core is as high as 2.5 x 1O 3 M pc -2 , which is about an order of magnitude higher than is 
typical for nearby normal galaxies (Sakamoto et al. 1999a). This value is however a lower 
limit, because the core is unresolved even in our high resolution observations (1.9" x 1.6"), 
The total flux within the radius of 2" (152 pc) is S C o = 1-2 x 10 2 Jy km s -1 ; the corresponding 
gas mass is M gas = 3.0 x 1O 8 M . 

3.5. Comparison with Previous Results 

On the molecular disk of NGC 3079, Sofue & Irwin (1992) identified three main compo- 
nents: an outer disk, spiral arms, and a central compact core. Recently Sofue et al. (2001) 
further resolved the central compact core, and classified it into two distinct components, 
a nuclear molecular disk and ultra-high-density molecular core, based on their appearance 
on maps and PV diagrams. In this paper, we re-classified the observed features into four 
distinct components, based on our new observations and analyses, and partly on theoretical 
considerations in §4. The outer disk and spiral arms in Sofue & Irwin (1992) correspond 
to our main disk and spiral arms respectively, and the ultra-high-density molecular core in 
Sofue et al. (2001) corresponds to our nuclear core. The nuclear molecular disk in Sofue et 
al. (2001) is different from our nuclear disk, and includes both our nuclear disk and main 
disk in part. 

We confirmed two peaks on a PV diagram in the central region, reported by Sofue 
& Irwin (1992), in Figure 10. This feature however cannot be confidently attributed to 
the nuclear ring (Sofue & Irwin 1992) because the two peaks may arise from our nuclear 
disk (see Figure 10 top and bottom) superposed on our nuclear core. Our central surface 
density E gas = 2.5 x 1O 3 M pc~ 2 is consistent with the value (2.2 x 10 4 M o pc~ 2 ) of Sofue 
& Irwin (1992) by considering the difference of the adopted conversion factors and applying 
the correction for inclination to the value of Sofue & Irwin (1992). We could not find 
unambiguous evidence for molecular outflow (Irwin & Sofue 1992) in our data. 

4. GAS DYNAMICS IN NGC 3079 

The near-perfect bisymmetric structures of the molecular disk are striking. These kinds 
of bisymmetry are often thought to be consequences of gas motions in a bar. Though a 
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stellar bar is not seen in optical/infrared photographs in this nearly edge-on galaxy, there 
is some indirect evidence for the presence of a bar: the peanut-shaped bulge (Combes et 
al. 1990; Shaw, Wilkinson & Carter 1993); the "figure-of-eight" pattern on the PV diagram 
(Kuijken & Merrifield 1995); the high central concentration of the molecular gas (Sakamoto 
et al. 1999b); the Ha velocity field well- fitted by an oval-orbit model in a weak bar (Veilleux, 
Bland-Hawthorn & Cecil 1999). We thus propose a model that NGC 3079 contains a weak 
bar. This model naturally explains the features in § 3: the main disk and spiral arms (or 
offset ridges) may result from gaseous rrx-orbits and associated crowding respectively. The 
nuclear disk may arise from gaseous rr 2 -orbits. This model however cannot interpret the 
nuclear core, which necessitates a massive component besides the bar at the galaxy center. 
Figure 11 shows schematic views of our interpretation, drawn on our model orbits (see a 
description in §4.1). We present below our detailed case for the existence of the weak bar, 
and discuss the origin of the main disk, spiral arms, nuclear disk and core. 



4.1. Damped Orbit Model 

We adopt a damped orbit model by Wada (1994) and Sakamoto, Barker & Scoville 
(2000) to describe gas motions in a weak bar potential. This model solves the equation of 
motion which includes a damping-force term to emulate the collisional nature of the gas, and 
obtains closed orbits for gas in a bar. This model has been applied to a galaxy NGC 5005 
(Sakamoto, Barker & Scoville 2000), and similar models have been successfully used to infer 
the presence of bars in edge-on galaxies which show characteristic patterns on PV diagrams 
(Binney et al. 1991; Achtermann & Lacy 1995; Garcia-Burillo & Guelin 1995; Kuijken & 
Merrifield 1995). 

Figure 12 displays our model of gas orbits in a bar potential that reproduces the main 
features in the zeroth-, first-moment maps, and PV diagram of NGC 3079. We here took 
the same potential model as Sakamoto, Barker & Scoville (2000) of 

$(r,0) = (1 -e cos 20)|log(l + (r/a) 2 ), (3) 

where the first and second terms stand for disk and bar potentials respectively, and 9 is 
the azimuthal angle from the bar major axis. For this particular case we adopted the 
characteristic radius a = 6" (456 pc), rotation velocity v = 240 km s -1 and bar strength 
e = 0.04 to match the observations. We chose the pattern speed of f2 b = 55kms _1 kpc _1 
based on consideration given in § 4.2. The inner and outer ILRs and corotation occur at the 
radii of 7.7" (582pc), 10.5" (798 pc) and 57.1" (4.34kpc) respectively. We assumed that the 
position angle of the bar (dashed line) from the north is P.A. int = 135° intrinsic to the disk, 
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and P.A. sky = 158° on the sky. On the galaxy disk, this position angle of the bar differs by 
30° from that of the galaxy major axis which runs horizontally in Figure 12 (top). 

Our parameters of the bar chosen to match the CO observations are in good agreement 
with the values that Veilleux, Bland-Hawthorn & Cecil (1999) have derived by fitting an 
oval-orbit model to their Ha velocity field of NGC 3079 (fib ~ 60kms _1 kpc -1 ; P.A. bar = 
130 ± 10°; axes ratio of oval orbits 0.7). Beyond our particular choices of the parameters, 
the features of gaseous closed orbits described below are general (Wada 1994; Lindblad & 
Lindblad 1994), and have been established in numerical simulations (Wada 1994; Wada & 
Koda 2001; Koda & Wada 2002). 

Closed orbits in a bar have an elongation in general; the position angle of the elongation 
if) is measured clockwise from the bar major axis (dashed line in Figure 12 top), in the rest of 
the discussion. For collisionless star-like particles, two families of closed orbits are dominant: 
rri-orbits have an elongation with ip = 0° and dominate outside of the outer ILR and inside 
of the inner ILR, while X2-orbits elongate with tp = 90° and are dominant between the 
ILRs. Orbits for gaseous particles show similar families at the same loci (Figure 12 top), 
while gaseous ^i-orbits have an elongation toward the leading edge of the bar (ip > 0°), and 
gaseous x 2 -orbits elongate opposite (ip < 90°), owing to the damping force. 

4.2. The Origin of The Main Disk and Spiral Arms 

The main disk and spiral arms may originate from gaseous rri-orbits. The position angle 
ip of orbits changes with radius owing to the damping force, producing crowded regions of the 
orbits (Figure 12). The crowding of the gaseous ^x-orbits which appear at the leading edge 
of the bar result in the spiral arms (or so-called offset ridges), which are similar to those on 
the molecular disk of NGC 3079. The slightly enhanced CO emission in the zeroth-moment 
map of NGC 3079 (Figure 2) may arise from the high gas densities in the crowded regions. 
The streamlines belonging to the a;i-orbits extend entirely to the disk, producing the main 
disk. 

Our model streamlines have sharp turns at the spiral arms. Such features, accompanied 
by shocks, have been observed by other authors (Athanassoula 1992; Athanassoula & Bureau 
1999) in full hydrodynamical simulations. These turns can produce double-peaked features in 
line-profiles observed along the spiral arms (Figure 9), and agree spatially with the S'-shaped 
twists on the first-moment map of NGC 3079 (Figure 2). Figure 13 shows a comparison 
between the observed and model PV diagrams; contours are from the observed one, while 
the symbols show the model one, where different symbols are used for gaseous particles on 
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the spiral arms (squares) and on the main disk (dots). On this diagram again, the main 
disk and spiral arms of NGC 3079 agree well with the gaseous rrx-orbits and their crowding 
respectively. This match between the observations and the model suggests that the main disk 
and spiral arms result from gaseous a;i-orbits and their associated crowding, respectively. 

In the context of our model, gas streamlines on the downstream side of the sharp 
turns are nearly perpendicular to our line-of-sight. Thus the line-of-sight velocity of the gas 
indicates the pattern speed of the bar. We estimate the pattern speed from the gradient of the 
spiral arms on the observed PV diagram, which appears to be = 55± 10 km s _1 kpc -1 ; this 
is our adopted value to describe gas orbits, and is consistent with the result of other authors 
(60kms _1 kpc -1 from Veilleux, Bland-Hawthorn & Cecil 1999, where the coincidence of 
corotation and the bar-end was assumed). 

4.3. The Origin of The Nuclear Disk 

The nuclear disk may result from gaseous X2-orbits. In Figure 13, the model PV diagram 
shows that the gaseous X2-orbits (crosses) are in good agreement with the nuclear disk of 
NGC 3079. Figure 14 displays PV diagrams along the galaxy major axis for the orbits 
within the outer ILR. The three panels show PV cuts with 10"-width (intrinsic to the galaxy) 
and positional offsets for the east, central and west sides of the model galaxy respectively. 
These panels reveal that the bends and twist of the ridges discussed in §3.3 (see Figure 10) 
are also reproduced in our model: the top and bottom panels show bends downward and 
upward respectively, while the middle panel shows a twist upward (at the left-hand side) and 
downward (right-hand). Since a purely rotating disk cannot produce this type of symmetric 
bends and twist, the gas on the nuclear disk must take oval orbits. Our gaseous a^-orbits 
naturally reproduce the observed properties of the nuclear disk. 

For later discussion in § 4.4, we mention the position angle of the a^-orbits on the 
nuclear disk here in advance. The major axis of the orbits are inclined counterclockwise 
to the south from our line-of-sight in our current configuration (Figure 11). On the other 
hand, orbits inclined in the opposite direction can also reproduce the nuclear disk on the PV 
diagram (Figure 4). These orbits however make the ridges on the eastern and western PV 
diagrams (Figure 14) bend in the opposite directions, and thus cannot match the observations 
(Figure 10). Therefore the orbits on the nuclear disk must have an inclination southward 
(counterclockwise) on the sky from our line-of-sight, like our gaseous X2-orbits. 

Many barred galaxies show an intense concentration of CO emission in their central 
parts (Sakamoto et al. 1999b; Sakamoto, Barker & Scoville 2000). NGC 3079 also has such 
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concentration in the nuclear disk region (see Figure 8). This feature is well explained in the 
scenario of gas motions in bar: gas flow on Xi-orbits like ours leads to offset shocks at the 
leading edge of a bar, and results in an inflow of gas toward the nuclear region. Then the 
inflowing gas on the Xi-orbits collides with gaseous rr 2 -orbits near the pericenter of the X\- 
orbits, entering onto the gaseous a^-orbits, thus increasing mass densities in the nuclear disk. 
This scenario has been confirmed in hydrodynamical simulations (Piner, Stone & Teuben 
1995; Athanassoula & Bureau 1999). 

4.4. The Origin of The Nuclear Core 

Our current model naturally explains the main features of the molecular disk of NGC 
3079, which are the consequences of gas motions in a weak bar. However the model does not 
explain the nuclear core which appears in the PV diagram as high velocity-peaks near the 
galaxy center (Figure 4). Owing to our ability to match the observations with our model for a 
bar, we attribute the high velocities to another origin, a massive central component, besides 
the bar. There still is a small possibility that the high velocities may arise from gaseous 
^i-orbits within the inner ILR which could elongate along our line-of-sight. However the 
gaseous rrx-orbits are not the likely origin for two reasons. 

(1) Our gaseous rri-orbits within the inner ILR (7.7") look near circular, thus cannot 
show high velocities at their pericenter. Numerical experiments have also shown similar 
features: circular orbits in the innermost regions (Piner, Stone & Teuben 1995; Athanassoula 
& Bureau 1999). In the intensive study of PV diagrams for galaxies with bars by Bureau & 
Athanassoula (1999) and Athanassoula & Bureau (1999), the Xi-orbits within the inner ILR 
do not produce any remarkable features (and velocity-peaks) on the PV diagrams. These all 
imply no significant contribution of the Xi-orbits within the inner ILR to our PV diagram. 
(2) Even if the gaseous rri-orbits are elongated, they must be near perpendicular to our 
line-of-sight, and cannot show us high velocities which occur at their pericenter, because 
the position angles ip of x±- and rr 2 -orbits, measured clockwise from the bar to the north 
(Figure 12 top), are intrinsically ip Xl = 0° and tp X2 = 90° respectively for collisionless stars. 
Owing to the damping force (working as viscosity), the two orbits drag each other and make 
■ip Xl > 0° and ip X2 < 90° for collisional gas. This mechanism does not allow ip Xl to become 
larger than ip X2 . In addition to ip Xl < ip X2 , our gaseous X2-orbits (nuclear disk) of NGC 3079 
have to be inclined southward (counterclockwise) on the sky from our line-of-sight for the 
reason discussed in § 4.3. Thus the elongated gaseous Xi-orbits (if they exist) must be near 
perpendicular to our line-of-sight, and cannot show the observed high- velocity features. 

There are many disk galaxies which exhibit double bars, i.e. global and nuclear bars (e.g. 
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Shaw et al. 1995; Friedli et al. 1996; Erwin & Sparke 1999). If NGC 3079 has a nuclear bar 
other than our global bar, it may result in the observed high-velocity features. This however 
is not the likely origin because such a nuclear bar is also thought to be a consequence 
of the X\- and rr 2 -orbits in our global bar, whose orientations cannot be aligned with our 
line-of-sight. Recently Heller, Shlosman & Englmaier (2001) discussed that a nuclear bar 
could be dynamically decoupled with a global bar in an evolutional phase, having a random 
orientation from the global bar. In that phase the x±- and x 2 -orbits have dissapeared (Heller, 
Shlosman & Englmaier 2001), which is not the case for NGC 3079. 

Central components with high velocities are often found in rot at ion- curve studies of 
external galaxies (Sofue et al. 1999). They are attributed to central massive components 
(Sofue et al. 1999; Takamiya & Sofue 2000) while gas streaming motions in a bar could 
be another possible origin (Sakamoto et al. 1999a). In this section, we clarified detailed 
kinematics and the effect of a bar on the gas disk of NGC 3079, and discussed that the high 
velocity should be attributed to a massive component at the galaxy center. 

5. CENTRAL ROTATION CURVE AND MASSIVE CORE 

We discussed in §4.4 that the high velocity peaks at the radius of ~ 1" in the PV diagram 
are not likely to result from noncircular motions driven by a bar, but are a consequence of 
circular rotation around the nuclear core. We here estimate the central dynamical mass 
using a rotation curve. 

5.1. Central Rotation Curve 

The observed PV diagram in the galaxy center suffers from large intrinsic velocity 
dispersions, and spill-in of emission from outer low velocity components owing to the finite 
spatial resolution and slit width. Thus the central rotation curve is not straightforwardly 
obtained from an observed PV diagram. We adopt the Takamiya & Sofue method (TS 
method; Takamiya & Sofue in private communication, and see Sofue & Rubin 2001) to obtain 
the rotation curve of the nuclear core. This method is designed to iteratively determine a 
rotation curve in an innermost region of a galaxy so that a model PV diagram calculated from 
the rotation curve and the galaxy's emission profile could reproduce the original diagram. 

Figure 15 displays PV diagrams and rotation curves of the nuclear core: the left-hand 
panel shows the observed one along the galaxy major axis with 1" slit-width, while the middle 
and right-hand panels show the calculated ones in the TS method and in the peak-tracing 
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method respectively. The emission profile at the center (~ 1") cannot be derived directly 
from our data owing to our spatial resolution (~ 2"). But the observations are indicative 
of its very sharp gradient, because even a point source can reproduce our observed profile 
through a convolution with the synthesized beam. Thus we assumed the emission profile 
as an analytic function of an exponential with the scale length of 1", which well reproduces 
observed features on the PV diagram. We assumed intrinsic velocity dispersions of 60kms _1 
(e.g. 15 — 50kms _1 in the Galactic center from Bally, Stark & Wilson 1988, which is less 
active than NGC 3079). 

The PV diagram and rotation curve from the TS method show that the envelope of the 
observed PV diagram requires a sharp rise of the rotation curve at the central region. The 
observed emission at around Okms -1 , which is not well reproduced in the TS method, may 
be contaminated by the emission from the nuclear disk surrounding the nuclear core. For a 
comparison, the right-hand panel displays the PV diagram using the rotation curve derived 
by tracing peak- intensity velocities of the observed PV diagram (the peak-tracing method); 
this method is often used to obtain rotation curves, but fails in reproducing the observed 
high velocity at the central part. 



Our spatial resolution is about 2" (152 pc), and the rotation velocity at this radius is 
220 km s" 1 from the rotation curve in the TS method (Figure 15 middle). Thus the dy- 
namical mass within the central 152 pc of NGC 3079 is estimated to be M dyn = 1.7 x 10 9 M Q . 
The ratio of dynamical mass to the total molecular gas mass (3.0 x 1O 8 M from §3.4) in the 
nuclear core is M gas /M dyn = 18%. 

The rotation curve is meaningful at the radius less than 2", even though the spatial 
resolution is only ~ 2": We can distinguish two emission peaks only if their spatial separation 
is larger than 2" on a single channel map or zeroth-moment map. However on the PV diagram 
and rotation curve, the additional parameter, velocity, provides a finer spatial resolution, 
because two emission peaks with different velocities can be resolved in the PV diagram even 
if their separation is less than 2". Then the spatial resolution will be determined by the 
typical error of a gaussian fit to a point source, which depends on the signal-to-noise ratio, 
and in our case is <C 1" (see Condon 1997). Therefore the velocity of 200 km s -1 at the radius 



5.2. 



Mass Estimation and Central Massive Core 



The dynamical mass within a radius r is estimated by 
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of 1" (76 pc) is meaningful, and implies a dynamical mass of M dyn = 7.0 x 10 8 M Q within the 
radius of 76 pc. The derived masses are listed in Table 4. 

Recently central supermassive black holes with mass 1O 7_8 M have been found in many 
AGNs (Miyoshi et al. 1995; Wandel, Peterson & Malkan 1999; Ishihara et al. 2001). In the 
case of NGC 3079, VLBI observations of H 2 maser emission indicate that its central mass 
within 0.5 pc is ~ 1O 6 M (Trotter et al. 1998). Thus our mass - 1O 9 M within - 100 pc 
cannot be attributed to this central supermassive black hole. The connection between this 
large dynamical mass and the nuclear outflow in parsec-scale (Henkel et al. 1984; Irwin & 
Seaquist 1988; Trotter et al. 1998) and kiloparsec-scale (Duric & Seaquist 1988; Ford et al. 
1986) in this galaxy will be an interesting issue for future investigations. 

6. CONCLUSIONS 

We have made CO (1-0) observations of the Ha/Radio lobe galaxy NGC 3079 with the 
Nobeyama Millimeter Array, and reported characteristic features and gas dynamics in the 
molecular disk in the central 4.5 kpc (1'). 

1. Our observations show four distinct components in the molecular disk: the main 
disk, spiral arms, nuclear disk, and nuclear core. 

2. The main disk extends along the galaxy major axis. We detected the inner region of 
this component, within the radius of ~ 2 kpc on a zeroth-moment map, while its full extent 
is not covered in our synthesis observations. Molecular gas shows a smooth distribution on 
the main disk, and has a gas mass of 5 x 1O 9 M within the central ~ 2 kpc. 

3. The spiral arms are superimposed on the main disk, forming an inverted S'-shaped 
pattern on the sky. Abrupt velocity changes of up to ~ 200 km s^ 1 are observed along the 
spiral arms, in S'-shaped twists of isovelocity contours and double velocity-peaked features 
on the spectra. This component and the nuclear disk form the so-called "figure-of-eight" 
pattern on a PV diagram, which sometimes appears in edge-on barred galaxies. 

4. The nuclear disk, appearing in PV diagrams, has an extent of ~ 600 pc radius. The 
molecular gas is intensely concentrated on the nuclear disk, where the gas surface density is 
an order of magnitude higher than that typical of nearby galaxies. Its appearance (bends 
and twist) on PV diagrams is indicative of oval motions in the gas, rather than circular. And 
the position angle of the oval orbits is inclined southward from our line-of-sight. 

5. The nuclear core has a radius smaller than ~ 150 pc, which is unresolved in our 
current observations. The gas mass within the central ~ 150 pc amounts to 3 x 1O 8 M . This 
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component shows a very high velocity ~ 200 km s 1 at the radius of the central ~ 100 pc on 
the PV diagram. 

6. A weak bar model explain the features of the main disk, spiral arms, and nuclear disk 
in zeroth- and first-moment maps, PV diagrams, and the spectrum at each spatial point. 
The main disk and spiral arms result from gaseous Xi-orbits and their crowding respectively. 
The nuclear disk arises from gaseous a^-orbits. The gas concentration in the nuclear disk 
is also explained in the context of our model: the gas on xi-orbits moves along spiral arms 
(or offset shocks), collides with the gas on X2-orbits, and is accumulated onto the nuclear 
disk. Assuming that the gas in spiral arms moves along streamlines perpendicular to our 
line-of-sight, the pattern speed of the bar can be estimated to be 55 ± lOkms -1 kpc -1 . 

7. Our model for a bar however does not explain the high velocity of the nuclear core. 
Moreover, any orbit caused by a bar is not likely to produce that component. Thus we at- 
tribute it to a central massive core with a dynamical mass of 10 9 M Q within the central 100 pc, 
which is three orders of magnitude more massive than the mass of a central supermassive 
black hole. 
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Table 1. PROPERTIES OF NGC 3079 



Parameter 


Value 


Reference 


Morphology 


SB(s)c sp 


1 




Sc pec: 


2 


Nuclear activity 


LINER 


3 




Type 2 Seyfert 


4 


Position of nucleus: 




5 


a(B1950) 


9 58 35.02 




5(B1950) 


55 55 15.4 




£>25 x d 2 5 


7.9 x 1.4 


1 


B T (mag) 


11.54 


1 


P.A. (°) 


165 


1 


Inclination (°) 


77 


6 


Distance (Mpc) 


15.6 


5 


Linear scale (pc arcsec -1 ) 


76 




-S'eO/im (Jy) 


50.17 ±0.054 


7 


-SloOjum (Jy) 


103.40 ±0.154 


7 


Lfir 


2.1 x 10 10 





a L Fm was calculated as 3.75 x 10 5 (£>/ Mpc) 2 (2.58S , 60Mm + 
Sioo^m) in units of L G . 



Note. - Reference.- (1) de Vaucouleurs et al. 1991; (2) 
Sandage & Tamman 1981; (3) Heckman 1980; (4) Ford et 
al. 1986; (5) Sofue & Irwin 1992; (6) Young et al. 1995; (7) 
Soifer et al. 1989 
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Table 2. PARAMETERS OF CO (1-0) CUBE 



Parameter Value 



Configuration AB+C+D 

Weighting NA 

Field of view (") 65 

Synthesized beam: 

major axis (") 1.9 

minor axis (") 1.6 

PA. (°) 105 

Velocity resolution ( km s _1 ) 10.4 

rms (mJybeam -1 ) 12 



Table 3. KINEMATIC PARAMETERS OF THE MAIN DISK 



Parameter Value 



Dynamical center: 

a(B1950) 9 58 35.00 

5(B1950) 55 55 15.90 

PA (deg) 169 

Inclination (deg) 79 

^(kms- 1 ) 1147 
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Table 4. GAS AND DYNAMICAL MASSES 



RADIUS a 










M dyn 




M gas /M dyn 


(arcsec) (pc) 


(1O 9 M ) 


(1O 8 M ) 


(%) 


26 2.0 x 10 3 


22 


46 (26) 


21 (12) 


2 1.5 x 10 2 


1.7 


3.0 (1.7) 


18 (10) 


1 76 


0.7 







a Radius in the galactic plane 

b Gas mass is calculated with a conversion factor 
X co = 1.8 x 10 20 cm- 2 [Kkms" 1 ]- 1 and the H abundance 
X = 0.707. The values in brackets are the ones with 
X co = 1.0 x 10 20 cm~ 2 [Kkms -1 ] -1 for high metallicity 
regions such as galactic centers (Arimoto, Sofue and Tsu- 
jimoto 1996) 
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Fig. 1. — Optical images of NGC 3079. (a) R-band image from the Digitalized Sky Survey, 
(b) HST WFPC2 images of J-band (grey scale) and Ha + [Nil] line emission (contour), 
obtained from the HST data archives (P.I. G. Cecil). Contours are at 2, 5, 15, 20, 30, 40, 
50, 60, 70, 80, 90, 100% of the Ha + [Nil] peak intensity. The circle represents the primary 
beam size of our CO (1-0) observations (65" HPBW). The absolute positions of both panels 
were calibrated using the USNO-A2.0 catalog (Zacharias et al. 2000), and uncertain to 0".5. 
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RISHT ASCENSION (B1BH) RIGHT ASCENSION (B19S0) 



Fig. 2. — CO (1-0) zeroth- and first-moment maps in the central 20" x 48" (1.5 kpc x 3.7 kpc) 
region of NGC 3079. The synthesized beam (1.9" x 1.6") is shown in the lower left corners. 
Left: The zeroth-moment map (integrated-intensity map). Contour levels are 1.5, 3, 5, 9, 
12, 15, 20, 25, 30, 40, 50, 60, 70, 80, 90, 100 % of the peak integrated-intensity 1.25 x 
10 2 Jybeam _1 kms" 1 . The clip level was set to 2.5a where la = 12mJybeam _1 . The 
primary beam attenuation has not been corrected. Right: The first-moment map (intensity- 
weighted mean velocity map). Contours are drawn with an interval of 50 km s" 1 ; The contour 
of 1150 kms -1 (V sys = 1147 kms" 1 ) is pointed to by an arrow. The clip level was set to 4a 
to make this map. 
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BIGHT ASCENSION Cfll 34-0) 

Fig. 3. — Channel maps in CO (1-0) emission. The channels have an interval of 10.4 km s -1 . 
Their central velocities (Vlsr in kms -1 ) are labeled at the lower-right corners. Crosses show 
the position of the dynamical center (Table 3). Contour levels are -3, 3, 5, 8, 12, 16, 20, 25, 
30<7, where la = 12 mJy beam -1 . Negative contours are dotted. No primary beam correction 
has been applied. 
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Fig. 3. — continued. 
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Fig. 3. — continued. 
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Fig. 4. — PV diagram of CO (1-0) emission along the major axis of NGC 3079 (P.A. = 165°). 
The slit-width was set to be 12", containing almost the entire emission in Figure 2. The 
axes are labeled relative to the dynamical center and systemic recession velocity (Table 3). 
Velocities have been corrected for inclination (77°). Contours are at 5, 10, 20, 30, 40, 50, 60, 
70, 80, 90, 100% of the peak intensity. 
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Vlsr [km* 1 ] 

Fig. 5. — Comparison of our CO (1-0) line spectrum with the single-dish result of IRAM 30 
m from Braine et al. (1993). The NMA cube was corrected for primary beam attenuation, 
convolved with the single-dish beam size (23") of IRAM, and sampled at the pointing centers 
of the IRAM observations (a w50 , 5 1950 ) = (09 h 58 m 35 s .4, +55°55 / ll".0). The gray line show 
the spectrum from the IRAM observations, whereas the black line show that from the NMA 
cube. 
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Fig. 6. — Schematic illustration of PV diagram. Four distinct components exist in the 
molecular disk of NGC 3079: main disk, spiral arms (or offset ridges), nuclear disk, and 
nuclear core. 
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Fig. 7. — Comparisons of CO maps (contours) with the Ha + [N II] map (grayscale; from 
the HST archives). The absolute position of the HST image was calibrated with the USNO- 
A2.0 catalog, and is accurate to about 0.5". The ellipses at the lower-left corners show the 
synthesized beam for the CO observations. The two panels show two sets of CO zeroth- 
moment maps with different integrated velocity ranges. Left: The CO map (contours) is 
derived by integrating the cube in the full velocity range where CO emission is detected. 
Contour levels are 3, 5, 9, 12, 15, 20, 25, 30, 40, 50, 60, 70, 80, 90, 100 % of the peak 
integrated intensity 1.25 x 10 2 Jybeam _1 kms -1 . The molecular disk coincides with the 
void of HII regions at the center, and its center is coincident with the root of the Ha lobe. 
Right: The CO map (contours) is derived by integrating the cube in the velocity range where 
the spiral arms are dominant in emission (AV = ±100 km s" 1 ). Contour levels are 5, 10, 20, 
30, 40, 50, 60, 70, 80, 90, 100 % of the peak integrated intensity 6.52 x 10 Jybeam" 1 kms" 1 . 
The spiral arms extend from the root of the Ha lobe, and wind counterclockwise. 
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Fig. 8. — Radial distribution of CO emission. CO integrated intensities are azimuthally 
averaged with corrections for inclination and primary beam attenuation. No correction for 
missing flux (about 10 — 30%) has been applied. The surface densities on the right-hand axis 
are calculated by assuming the Galactic conversion factor Xqo = 1-8 x 10 20 cm" 2 [K km s -1 ] -1 
(Dame, Hartmann and Thaddeus 2001). The two horizontal bars show the beam size along 
the major and minor axes (1".9 x 1".6). The solid line presents the derived emission profile, 
while the dotted line shows a fitted exponential profile, i.e. 5.97exp(— r/7".24), at radii of 
r > 5". 
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Fig. 9. — CO (1-0) line profiles and zeroth- moment map for reference. The line profiles are 
sampled at 5 x 7 points on the 3"-spacing grid (crosses on the left-hand panel), centered at the 
dynamical center (Table 3). Double-peaked features in the line profiles are observed along 
the spiral arms: arrows point to two examples. The double-peaks originate from distinct 
velocities between the main disk and spiral arms (or offset ridges). 
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Fig. 10. — PV diagrams along the major axis (P.A. = 165°) with offsets along the minor 
axis of +2", 0" and —2" ("+" denotes east). The slit-widths are 1". The axes are labeled 
relative to the dynamical center and systemic recession velocity (Table 3). Velocities have 
been corrected for inclination (77°). Contours in all the panels are at 2, 3, 4, 5, 6, 7, 8, 9, 
10 times 4.0 x 10 -2 Jy beam -1 . 
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Fig. 11. — Schematic views of our interpretation of the molecular disk in NGC 3079. Four 
distinct components, the main disk, spiral arms/offset ridges, nuclear disk, and nuclear core, 
are drawn on the streamlines of our model orbits in the bar reference frame (§4.1). The 
bar has the position angle of 135° intrinsic to the galaxy from the north. Left: Face-on 
view. The galaxy major axis lies horizontally, while the bar runs along dashed line which is 
inclined clockwise by 30° from the horizontal. Right: Projection on the sky for P. A. = 165° 
and i = 77°. 




Fig. 12. — Model gas motions in the molecular disk of NGC 3079. Top: Face-on view 
of model orbits in a weak bar potential, on the frame rotating with the bar. The galaxy 
major axis runs horizontally. The bar is rotating clockwise, and runs along the dashed line 
(inclined from the horizontal by 30°). The orbits are calculated using a damped orbit model 
(Wada 1994). The inner and outer ILR and corotation occur at the radii of 8", 11" and 57" 
(1" = 76 pc) respectively. Bottom: Model PV diagram, which cuts the model galaxy along 
the galaxy major axis, with the width containing the full region at the top panel. The solid 
line shows the circular rotation curve. Different symbols are used for different orbits: circles 
are xi-orbits inside the inner ILR, crosses are X2-orbits, squares are xi-orbits outside the 
outer ILR in front of the sharp turns on the streamlines, and dots are also the Xi-orbits at 
the rests. 
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Fig. 13. — Comparison of the model and observed PV diagrams. Contours trace the observed 
PV diagram (Figure 4) at 8 and 40% of the peak intensity. Symbols plotted are the same 
as those in the bottom panel of Figure 12. Our model reproduces the main features of the 
molecular disk of NGC 3079 except the nuclear core. This indicates the presence of a massive 
core. 




Fig. 14. — Model PV diagrams for the orbits within the outer ILR. The diagrams are cut 
along the major axis for east (y > 0"), center and west (y < 0") sides, where y is the vertical 
axis on Figure 12 (top). The cut has a 10" width in the galaxy's face-on plane. The symbols 
are the same as those in Figure 12 (bottom): crosses are a^-orbits, and circles are Xi-orbits 
inside the inner ILR. The model well represents the bends (east and west) and twist (center) 
discussed in §3.3 (see Figure 10). 
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Fig. 15. — PV diagrams and rotation curves. The slit is 1" wide. Contours are 20, 40, 60, 
80% of the peak intensity in each panel. Left: Observed PV diagram in the central region 
of NGC 3079. Middle: PV diagram and rotation curve derived in the Takamiya & Sofue 
method. Right: PV diagram and rotation curve derived in the peak-tracing method. 



